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SECTION  1 


INTRODUCTION 


The  governing  equations  in  fluid  mechanics  are  usually 
expressed  in  terms  of  Cartesian  coordinates.  They  are  not 
readily  convertible  into  forms  useful  for  practical  appli¬ 
cations  to  one-,  two-  and  three-dimensional  cylindrical  and 
spherical  problems.  One  of  the  purposes  of  this  report  is  to 
transform  the  mass,  momentum  and  energy  conservation  equations 
for  a  nonsteady  flow  of  a  compressible  viscous  single-phase 
fluid  from  expressions  in  vector  and  dyadic  notations  to  those 
in  terms  of  orthogonal  cuvi linear  coordinates  and  then  to 
specialize  the  equations  for  Cartesian,  cylindrical  and 
spherical  coordinates. 

Another  purpose  of  this  report  is  to  provide  a  theoretical 
basis  for  usage  of  Eulerlan  sliding  grids  in  modern  computa¬ 
tional  fluid  mechanics  by  establishing  a  kinematic  transport 
theorem  which  is  a  generalization  of  the  Reynolds  transport 
theorem.  Based  upon  the  theorem  the  conservation  equations 
in  integral  form  are  first  formulated  in  terms  of  moving 
volume  and  surface  elements  and  then  converted  into  the  corres¬ 
ponding  equations  in  differential  form  by  using  the  general¬ 
ized  definition  of  the  vector  operator  V. 

In  order  to  present  an  overall  survey  of  the  basic  equations 
this  rojyrt  also  summarizes  the  equations  for  the  vorticity, 
entropy  enthalpy  and  Bernoulli  equation. 
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SECTION  2 


INTEGRAL  FORMULATION 

2.1  KINEMATIC  TRANSPORT  THEOREM. 

Let  x  =  (x1#  x2,  x3)  denote  the  rectangular  spatial  ("Eulerian") 
coordinates  which  identify  a  fixed  point  in  space.  Let  X  = 

(X^  X2,  X3)  denote  the  rectangular  material  ("Lagrangian“) 
coordinates  which  identify  a  fluid  particle  in  motion.  Let 
F  =  F  (x,  t)  represent  any  arbitrary  single-valued  scalar  or 
vector  point  function  (of  position  x  and  time  t)  possessing 
continuous  derivatives.  The  function 

M^  =  J  F(x  ,  t)dV  «  Jv  F [x  (X,  t),t]J  dV^  ,  (1) 

V  ° 

where  V  =  \7(t)  denotes  a  material  volume  (that  is*  a  volume 
moving  with  the  fluid) ,  is  a  well-defined  function  of  time. 

In  Eq.  (1)  the  Jacobian  (Ref.  1,  p.33) 

dt7  aCXwX-jX-,)  /Sx.v 

J  -  xra  ■  35 "  det(s5g<  (i  *  1-2-3‘  *  ■ l'2'31  (2> 

relates  the  element  dv  of  the  moving  volume  V  in  the  x-variables 
to  the  element  dVQ  of  the  fixed  volume  VQ  »  \f(o)  *  v(t)  as  t  *  0 
in  the  X- variables. 

Using  Euler ' s  expansion  formula  (Ref.  2) 

=  JV • v  (3) 
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(where  v  is  the  velocity  vector  of  the  fluid  motion)  and  the 
relation  between  the  material  time  derivative  and  the  spatial 
derivatives  (Ref.  3,  Eq.  (3.6)) 


aS  *  Jt  +  Z’7  (4> 

we  can  express  the  material  time  derivative  of  as 

H)dvo 

o 

"/  [H  +  P(7’^  Jdvo 
Vo  1 

“/I  [H  *  v*VF  ♦  F(V*v) 

[|f  +  aiv  (p  v>  d?  .  cs> 

Then  in  view  of  Green's  transformation  (Ref.  1,  Eq.  (7.2)) 
for  any  vector  or  tensor  field  $ 

f  div  $  dV  =  f  dA  «  $  ,  (dA  *  ndA)  (6) 

we  obtain 
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(7) 


/v  ^ 

or 

Jt  /yFdV  *  "  Fv*ndA  +  ,  (7a) 


where  V  denotes  the  volume  fixed  in  space  which  instantaneously 
coincides  with  the  material  volume  V ,  A  denotes  the  surface 
bounding  the  volume  V  and  n  denotes  the  unit  vector  along  the 
outward  normal  to  A.  Eq.  (7a)  is  known  as  the  transport 
theorem  of  Reynolds  (Ref.  lf  Eq.  (25.4)  and  Ref.  4,  §  14). 

It  is  a  kinematic  relation  independent  of  any  meaning  attached 
to  F.  All  fluid  physics  is  contained  in  the  term* 

Now  let  us  consider  a  volume  V  «  V(t)  with  bounding  surface 
A  «  A(t)  sliding  with  respect  to  the  fluid  and  obeying  the 
relation 

<dV)  «  unidA  ,  (8) 

where  u  is  the  local  surface  velocity.  We  can  write 

k  ■  fjf a-aa*  .  »> 

V  A 

To  each  time  instant  t  when  the  fixed  volume  v  coincides  with 
the  material  volume  V  there  corresponds  a  time  instant  t.  when 

I-  A 

the  moving  volume  V  coincides  with  the  same  material  volume 
$  Instantaneously.  At  t^  Eq.  (d)  becomes 

jfpj  ^dV)  -jTPiUi.ndA  ,  (10) 
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where  =  F(x,t1),  u^  =  n(x,t^)  and  (dV)  =  -^(dV) 

Since  V,  A  and  n  in  Eq.  (10)  have  the  same  meaning  as  those 
in  Eq.  (7)  and  dV  and  dA  (at  t)  in  Eq.  (7)  equal  to  d^  and 
dk  (at  t^)  in  Eq.  (10)  we  may  add  Eqs.  (7)  and  (10)  to  obtain 

+  F1 il^>]  -  -  &  <n  -  Fi-i*  -s»  +  ^  «v  •  <ll> 


In  case  V  and  V  can  coincide  with  V  at  the  same  time 
(Cor  example,  if  =  0  or  *  v)  Eq.  (11)  becomes 

^(Pdv)  -  -  ^Ftv  -  u)-ndA  +  ^  , 

or  (since  V  is  fixed  in  space) 

jg  J FdV  «  -  <J^F(v  “  *«dA  «•  ,  (Ub) 

This  is  a  generalization  of  the  Reynolds  transport  theorem 
and  gives  the  fundamental  transport  formula  for  the  time  vari¬ 
ation  of  the  volume  integral  of  any  fluid  state  quantity  F 

over  a  moving  volume  which  instantaneously  coincides  with 

n, 

the  material  volume.  As  u  *  0  (V  =  V,  hence  «  t)  Eq.  (lib) 

reduces  to  Eq.  (7a)  which  is  the  basis  of  the  Eulerian  grid 

a  _ 

representation  in  numerical  algorithm.  As  u  *  v  (V  =s  V,  hence 
t^  =  t)  Eq.  (lib)  becomes  a  statement  of  the  Lagrangian  mesh 
representation.  Ref.  5  shows  an  example  of  the  application 
of  Eq.  (lib). 

2.2  CONSERVATION  LAWS. 

The  mass,  momentum  and  energy  conservation  equations  for  a  non¬ 
steady  flow  of  a  compressible,  viscous,  heat-conducting  fluid 


(tj^  *  t) 

(11a) 
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derived  and  discussed  in  reference  books  on  gas  dynamics  (see, 
for  example,  Refs.  6  to  13)  have  been  summarized  in  recent 
books  on  computational  fluid  mechanics  (see,  for  example. 

Refs.  14  to  16)  in  compact  vector  and  dyadic  (second  order 
tensor)  notations.  In  a  fixed  Eulerian  frame  of  reference 
these  equations  in  integral  form  (which  is  more  basic  than 
the  differential  form  in  expressing  the  physical  laws)  are  as 
follows : 


mass: 


n Xpdv  +  4  p--aa  =  °  ■ 


momentum:  J pvdV  +  ^pv(v«n)  dA  *  ^a«ndA  +j£fdV  > 

energy :  J* pEdV  +  J^pEv*ndA  =  ^  (a*v  -  c[)  *ndA 


t-j^f  *vdV  +  J 


GdV  . 


In  these  equations,  p  is  the  density,  E  the  total  specific 
energy  (''specific'*  means  per  unit  mass) 

E  =  e  +  ^  v^  (15) 


(where  e  is  the  specific  internal  energy) ,  a  is  the  stress 
tensor,  the  heat-flux  vector,  f  the  external  force  vector 
per  unit  volume ,  and  G  is  the  energy  generation  per  unit 
volume,  and  V  is  a  control  volume  fixed  in  space,  A  is  the 
surface  hounding  V,  and  n  is  the  unit  vector  along  the  out¬ 
ward  normal  to  A. 


In  these  equations  the  properties  of  the  fluid  need  not  be 
continuous  functions  of  space  and  time. 


By  comparison  of  Eq.  (7a)  with  Eqs.  (12)  to  (14)  the  function 
F  may  be  identified  as  p,  pv  and  pE  respectively  and  the  term 
dM^/dt  as  the  right-hand  side  terms  of  each  of  these  equations. 
Then  in  terms  of  the  moving  elements  (d^  =  dV  and  dA  =  dA 
instantaneously)  the  mass,  momentum  and  energy  conservation 
equations  in  case  t^  «  t  (exactly  or  approximately)  can  be 
expressed  exactly  or  approximately  as  a  single  vector  equation 
according  to  Eq.  (lib) : 


Eq.  (16)  provides  a  theoretical  basis  for  usage  of  Euler ian 
sliding  grids  in  finite-difference  numerical  algorithm  (see, 
for  example.  Ref.  5) . 

2.3  CONSTITUTIVE  RELATIONS. 

The  basic  dependent  variables  in  Eqs.  (12)  to  (14)  or  Eqs.  (16) 
are  p,  v  and  E  (or  e) .  Constitutive  relations  for  g  and  <j 
must  be  added  to  these  equations  in  order  to  obtain  a  closed 
system.  We  are  concerned  here  with  the  case  of  Newtonian 
fluids,  i.e.,  by  definition,  fluids  such  that  the  stress  tensor 


7 


is  a  linear  function  of  the  velocity  gradient.  From  this 
definition,  excluding  the  existence  of  districted  force 
couples,  results  Newton's  law,  also  called  the  Navier-Stokes 
law,  for  g: 

i  “  -  P|  +  l  ,  (18) 

with 

J  *  X  (7»v)I  +  2yD  (19) 

and 

(vv  +  <VZ)  *]  ,  (20) 

the  superscript  t  denoting  the  transpose  of  a  tensor.  In 
Eqs.  (13)  to  (20)  p  is  the  pressure,  j  is  the  viscous  (or 
deviatoric)  stress  tensor,  J  is  the  unit  tensor,  u  and  X  are 
the  first  (shear)  and  second  (dilatational)  coefficients  of 
viscosity,  and  £  la  the  tensor  of  rates  of  deformation. 
Furthermore,  the  fluid  is  assumed  to  obey  Fourier's  law  of 
heat  conduction  for  <j: 

a  *  -kVT  ,  (21) 

where  T  is  the  absolute  temperature,  and  k  is  the  thermal  con¬ 
ductivity  coefficient.  Many  fluids,  in  particular  air  and 
water,  follow  Newton's  law  and  Fourier's  law. 

2.4  TEEHMODVNAMIC/TBANSPORT  PROPERTIES. 

The  state  variables  p,  e,  p,  T  and  the  specific  entropy  S  are 
connected  by  thermodynamic  relations  (assuming  local  thermo¬ 
dynamic  equilibrium) . 


We  consider  the  case  of  a  simple  fluid  such  that  all  its  thermo¬ 
dynamic  properties  can  be  deduced  from  a  single  fundamental 
relationship  which,  for  a  compressible  fluid,  can  be  chosen  of 
the  type 


S  =  S (p,e)  .  (22) 

From  thi3  relationship  p  and  T  are  obtained  in  terms  of  the 
basic  variables  p  and  e  from 

P  ■  -»2l(lf)e  '  T  -  WSiT,  •  (23> 

An  important  special  case  is  a  perfect  gas  with  constant  spe¬ 
cific  heats  c^  and  ov.  For  ^uch  a  gas  the  laws  of  state  are 

p  ■  (Y-l)oe,  |  Y  =  cp/cv)  (24) 

and 

e  -  cvT  .  (25) 

The  viscosity  and  thermal  conductivity  coefficients  depend 
on  the  local  thermodynamic  state;  in  most  conditions  they  depend 
only  on  the  temperature: 

y  ■  U(T),  X  ®  A  (T)  ,  k  «  K,T)  •  (26) 


The  coefficient 


tc  *  3X  +  2p  (27) 

is  called  bulk  viscosity  coefficient  In  the  "Stokes  relation" 
(Ref.  3,  p.  238} 


3X  +  2p  “  0 


(2C) 


o 


it  is  assumed  to  be  zero.  However,  except  for  very  special 
conditions,  for  example,  monatomic  gases,  there  is  no  reason 
to  assume  3X  =  -2y.  (Ref.  9,  p.337  and  Ref.  7,  p.540.) 


-a  *c 


■a..  *1. 


SECTION  3 


CONSERVATION  LAWS  IN  DIFFERENTIAL  FORM 

3.1  CONVERSION  BY  GENERALIZED  DEFINITION  OF  OPERATOR  7. 

If  the  properties  of  the  fluid  are  continuous  and  sufficiently 
differentiable  in  some  domain  of  space  and  time,  then  the  con¬ 
servation  equations  in  integral  form  can  be  converted  into  an 
equivalent  set  of  partial  differential  equations  through  a 
generalized  definition  of  the  vector  operator  V  (Ref.  17,  p.40) 

74>  =  lim  ^  §.  4>ndA  ,  (29) 

v+o  v  “ 

where  4>  is  an  unspecified  (scalar,  vector  or  dyadic)  function 
of  position,  V  is  the  volume  enclosed  by  a  surface  A  to  which 
the  point  P  at  which  74)  is  to  be  calculated  remains  interior, 
while  the  largest  dimension  of  A  tends  to  zero,  and  the  multi¬ 
plication  in  V4>  may  be  scalar,  vector  or  dyadic.  In  particular 
when  4>  »  b,  a  vector,  the  scalar  product  becomes 

7«b  «  lim  hi.  b*ndA  ,  (30) 

“  V+0  V/A“‘ 

which  is  a  scalar ,  and  when  $  »  | ,  a  dyadic,  the  dot  product 
becomes 


7«|l  «  lim  |  |*ndA  ,  (31) 

which  is  a  vector.  Eqs.  (30) and  (31)  can  be  also  obtained  by 
dividing  Eq.  (6)  by  V  and  then  taking  the  limits  of  both  sides 
as  V  approaches  to  zero. 
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3.2  CONSERVATIVE  FORM  . 


Dividing  Eqs.  (12)  to  (14)  by  V,  taking  the  limit  of  every 
term  as  V  approaches  zero  and  applying  Eqs.  (30)  and  (31)  we 
obtain  the  conservation  equations  in  differential  form: 

mass:  +  7»(pv)  *  0  ,  (32) 

momentum:  (pv)  +  V • (pvv  -  g)  =  f  ,  (33) 

energy:  ^  (pE)  +  V»(pEv  -  g*v  +  <j)  ■  f«v  +  G  .  (34) 

An  alternate  form  of  the  energy  equation  Eq.  (34)  is  in  terms 
of  the  enthalpy  per  unit  mass 

h  •  e  +  p/p  (35) 

Substituting 

E  ■  h  +  j  V2  -  p/p  *  H  •  p/p  (36) 

into  Eq.  (34)  we  obtain 

(pH)  ♦  V*  (pHv  -  p v  -  g*v  <j)  ■  f*v  +  +  o  »  (37) 

where 


H  *  h  +  J  v2  a  E  +  p/p  (38) 

is  the  total  enthalpy  per  unit  mass  (which  is  also  known  as  the 
“specific  stagnation  enthalpy"  or  for  nonviscous  and  nonheat¬ 
conducting  fluid  in  steady  motion  the  "total  energy  per  unit 
mass" ,  f .  18,  p.  33) . 
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Substituting  Eg.  (36)  into  Eq.  (14)  we  obtain  the  integral 
form  of  the  conservation  equation  for  pH  in  terms  of  fixed 
elements  dV  and  dA 


TC XPHdv  +  /a  pH-- 


ndA 


1  (= •  !£~SL+P2^  *  2!^  +  fvll-Z+G*  HW 


(37a) 


The  corresponding  conservation  equation  in  terms  of  moving 
elements  dty  and  dk  in  case  t^  *  t  is  Eq.  (16)  with 


W  *  (pH) ,  I  *  (g.v-£+pv) ,  Q  *  (f  .V-H3+  |£)  (17a) 

3.3  NONCONSERVATIVE  FORM. 


Sqs*  (32)  to  (34)  are  mass,  momentum  and  energy  equations  in 
"conservative"  or  "divergence"  form.  (See  Section  III-A-3  of 
Ref.  14  for  the  meaning  and  beneficial  effects  of  using  this 
form) .  The  corresponding  equations  in  nonconservative  form 
in  which  g  follows  Newton’s  law  are  (Ref.  16,  Section  1.1): 

mass:  3t  ^  pV*v«’0  ,  (39) 

momentum:  •  i  ♦  *  -  *  +  u72v  +  (X+p)7(7*v) 

+  (7 *v)  7X  +  2g*7y  ,  (40) 

energy:  p  ^  +  p7»v  »  4-7*c[  +  G  ,  (41) 
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where  <$  is  the  dissipation  function 


«-  tjVv  *  A(V*v)2  +  2yD:D  (42) 

33  w  SB  3C 

and  ^  is  the  material  derivative  given  in  Eg.  (4). 

An  alternate  form  of  the  energy  equation  Eg.  (41)  is  in  terms 
of  the  specific  entropy  S  such  that  (Ref.  3,  Eq.  (33.1)) 

TdS  -  de  +  pd(|)  .  (43) 

By  use  of  Eqs.  (39)  and  (43)  we  obtain  from  Eq.  (41) 

pT  g|  -  « -  7«q  +  G  .  (44) 

Adding  Eq.  (32) (multiplied  by  S)  and  Eq.  (44)  (divided  by  T) 
we  have  the  conservative  form  of  the  conservation  equation 
for  pS 


(pS)  +  V.(pSv)  *  !  (*~V.a-K5)  .  (44a) 

Taking  the  integral  of  Eq.  (44a)  over  a  fixed  control  volume  V 
and  applying  Green's  transformation  Eq.  (6)  we  obtain  the 
integral  form  of  the  conservation  equation  for  pS  in  terms  of 
fixed  elements  dV  and  dA 

n  X  pSdV  +  y  pSv»ndA 


®jfy  (*-V.q+G)dV  (44b) 
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The  corresponding  conservation  equation  in  terms  of  moving 
elements  d$  and  d&  in  case  t^  =  t  is  Eq.  (16)  with 

W  *  (pS) ,  S  -  (0) ,  Q  *  ($/T-7 .g/T+G/T)  (17b) 

3.4  VORTICITY  DISTRIBUTION. 


The  vorticity  vector 

£  -  7  x  v  (45) 

gives  the  intrinsic  rotation  of  each  fluid  element. 

By  taking  the  curl  of  both  sides  of  Eq.  (40)  (divided  by  p) 
and  noting  that 


7  •  n  »  0 


(46) 


and 


7*7  (scalar  function)  *  0 


(47) 


we  obtain  the  general  equation  of  vorticity  distribution 
||  -  (Q*7)v  +  Q  (7»v)  +  7(i)  x  7p 


I  7  x  f  +  V(h  x  f  +  H.  72t)  -  V(£)  x  7  x  ft 
P  —  P  —  P  —  P  — 


+  7(~lii)x  7(7. v)  -  i  7X  x  7(7. v) 

+  (7.v)7(i)  x  7X  +  |  7  x  (g.7y)  +  7  (|)  x  (g.7y) 


(48) 


If  u  and  X  are  not  functions  of  space  variables  Eq.  (48) 
reduces  to 
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(49) 


+  (X+2ii)V(i)  *  7(7**)  -  U7(g)  xvxfl 

If  the  flow  is  incompressible  which  is  characterized  by  the 
condition 


7*v  »  0  (50) 

and  implies  that  p  is  constant  along  a  fluid  particle  trajec¬ 
tory  but  not  necessarily  independent  of  space  variables  (as 
in  stratified  flows,  Ref.  1$,  Section  1.3)  Eg.  (48)  reduces  to 

'  -  (£’?>v  +  7(|)  *x  7p 

•  |  7  x  f  +  V(i)  + 

-7(|)*7*n  +  |  7  x  (g.7u>  ♦  7(|)  x  (g.Vu)  .  (51)  , 

If,  in  addition,  p  is  constant  everywhere  Eq.  (51)  reduces  to 

dfi  tn  _ 

g£  -  (£*7)V 


f  +  £  v2n 

—  0  — 


X  7x0  +  £  7  X  (*5»7y) 


(52) 
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i |L  -.j •*-«  ,»>  V  V  Vr.rt’V *>:*  ■.*’*/ WW 


If,  in  addition,  y  is  not  a  function  of  space  variables 
Eq.  (52)  reduces  to  (Ref.  16,  Eq  .  (1.37)) 

g§  “  (0«7)v  =  i  7.x  f  +  H  72a  .  (53) 


This  equation  is  usually  associated  with  an  equation  for  a 
solenoid  stream- function  vector  ¥  such  that 

v  =  7  x  l  ,  (54) 

which  automatically  satisfies  the  incompressibility  condition 
Eq.  (50) .  The  equation  satisfied  by  V  is  derived  by  applying 
the  curl  operator  to  Eq.  (54)  and  using  definition  Eq.  (45) 
to  obtain 

+  ft  -  0  (55) 

since 

-  0  (56) 

as  V  is  solenoid. 

Let  us  return  now  to  consider  the  case  of  compressible  flow. 

For  an  ideal  fluid  (y*X*0)  either  Eq.  (48)  or  Eq.  (49)  reduces 

to  [Ref.  18,  Eq.  (7-20)] 

-  (ft*V)v  +  ft  (7«V) 

a  I  v  k  f  +  V(|)  5*  f  -  y  »  (i  7p)  .  (57) 

r  r*  r 
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Another  form  of  this  equation  can  be  obtained  by  using 
Eq.  (43)  [Ref.  18,  Eq.  (7-21)j  : 


-  (n»7)v  4-  fl(Vv) 

*  i  ?  x  f  4*  7(i)  )t  f  +  7T  x  7S  .  (58) 

If  p  is  a  function  of  p  only  the  last  term  of  Eq.  (57)  can 
be  written  as 


7  x  (i  7p)  -  7  x  ,  (59) 

which  is  zero  according  to  Eq.  (47).  If  the  specific  entropy 
S  is  a  constant,  then  the  last  tens  of  Eq.  (58)  is  zero 
because  7S  -  &.  Both  conditions  are  satisfied  if  the  fluid  is 
isentropic,  i.e.,  the  fluid  has  the  sane  entropy  everywhere. 

In  the  absence  of  any  external  force  f ,  if  ft  is  zero  for  one 
instant  at  every  point  of  the  flow  field  then  dfi/dt  - 
Therefore  the  vorticity  at  every  point  of  the  field  will 
remain  zero  and  the  flow  is  irrotational.  On  the  other  hand, 
if  the  flow  is  not  isentropic,  i.e.,  the  fluid  has  different 
entropy  at  different  points  of  the  field,  the  last  term  of 
Eq.  (56)  will  cause  the  vorticity  to  be  different  from  zero  in 
the  next  instant.  Therefore,  nonisentropy  flows  cannot  be 
irrotational.  Bence  irrotational ity  implies  isentropy,  but 
iaentropy  does  not  imply  irrotationality. 

Multiplying  Eq.  (48)  by  p,  expanding  the  term  (pn»7)v  by 
formula  (X)  on  p.  44  of  Ref.  17  and  adding  Eq.  (32)  (multiplied 
by  0)  we  have  the  conservative  form  of  the  conservation  equa¬ 
tion  for  pfi 
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(48a) 


(pn)  +  v» (pv  q)  -  -y*(pa  v)  +  Q± 


where 


£l  *  -v[7*  (pn)]  -  pfl(7»v)  -  p7(|)x7p  +  p7(i)xf  +  7xf 


+  U72ft  +  p7(^i-)x7(7‘v)  -  P7(H.)x7xn 


+  p(7*v)?(|)x7X  -  7Xx7(7*v)  +  27x(g*7y) 


♦  pV(^)x(g.7u) 


(48b) 


Taking  the  integral  of  Eg.  (48a)  over  a  fixed  control  volume  V 
and  applying  Green's  transformation  Eq.  (6)  and  Egs.  (6)  and 
(?)  on  p.  53  of  Ref.  1?  we  obtain  the  integral  form  of  the 
conservation  equation  for  pH  in  terms  of  fixed  elements  dV 
and  dA 


£  pftdV  +  f pfi(v*n)dA 
-/ pv(a*n)dA>  j^QjdV 


(48c) 


The  corresponding  conservation  equation  in  terms  of  moving 
elements  d$  and  dA  in  case  *  t  is  £c.  (16)  with 


w«(pn),  E  »  (pv  ft ) ,  Q  *  (Qj_) 


(17c) 


Combining  Eqs.  (17a) ,  (17b)  and  (17c)  with  Eqs.  (17)  we  con¬ 
vert  Eq.  (16)  to  a  vector  equation  representing  six  conser¬ 
vation  equations  in  integral  form  with 
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0 


0 


2* v  -  cr 

a*v  -  g.  +  pv 
0 


\ 


pv  n 


Q  = 


f 

f*v  +  G 


.  (17d) 


f*v  +  G  + 

-  ot 

$/T  -  7*g/T  +  G/Tj 


\ 


3.5  BERNOULLI  EQUATION. 


As  discussed  in  Ref.  18,  Section  A, 6,  aside  from  the  boundary 
layer,  the  effects  of  viscosity  and  heat  conduction  can  be 
neglected  for  the  majority  of  gas  dynamics  problems.  Further¬ 
more,  the  heat  addition,  except  in  problems  involving  combus¬ 
tion,  is  either  zero  or  ^ery  small.  Then  under  ordinary 
conditions  the  gas  behaves  very  much  like  an  ideal  fluid 
(i.e.,  a  nonviscous  and  nonhead-conducting  fluid).  Therefore, 
one  of  the  fundamental  problems  of  gas  dynamics  is  to  study 
the  adiabatic  flow  of  an  ideal  gas. 


Let  us  further  assume  that  the  density  of  the  fluid  is  a  func¬ 
tion  of  pressure  only  (i.e.,  the  fluid  is  barotropic,  Ref.  3, 
p.  150)  and  that  the  external  force  per  unit  mass  is  repre¬ 
sentable  by  a  potential^ 


f/p  =  -  Vtf  .  (60) 

Then  by  using  the  relation  (Ref.  17,  p.  44) 

(v*V)v  a  ^  Vv^  ~  v  X  (V  X  v)  (61) 

and  definition  of  £  (Eq.  (45))  Eq.  (40)  can  be  written  as 


'  ->  ■  •  „ »  •  ■  >  ■  V  ',‘V  V  .  -v. 


.  '.<  »< .  * 


V-V 


.•  v; 


v 
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-  V  x  a  =  —■  Vp  -  v(j  v2J  -  v<p*  -VX  ,  (62) 

where 

x  mf  i?  *  i  v'2  +  *  *  (63) 

Let  s  be  a  parameter  along  an  arbitrary  curve  in  space.  At 
any  point  (with  position  vector  r)  along  this  curve  dr/ds  is 
a  unit  tangent  vector  to  the  curve.  The  component  of  Eq.  (i2) 
in  the  direction  of  dr/ds  is 


3v  dr  dr  dr 

It  *  3s  *  *  3s  “  “  at  *  7X  * 


Since 


*  f{ 3v  dr\  3v  dr 

«  '  V  US  ‘  atj43  *  st  •  as 


we  may  express  Eq.  (64)  as  (Ref. 19,  p.  121) 

dr  dr 

(v  x  g)  *  gf  “  af  •  » 


where 


x.y^  + 1  v2  +  *+/( If  •  H)  as  . 


At  any  point  in  time  and  at  every  point  in  the  flow  field 
let  dr/ds  represent  a  tangent  vector  to  a  streamline  (in  the 
direction  of  v)  or  a  vortex  line  (in  the  direction  of  fi) . 


21 


Then  the  scalar  triple  product  on  the  left  side  of  Eq.  (66) 
is  zero  and  the  Bernoulli's  equation 

X  *  a  constant  (68) 

is  valid  along  any  streamline  or  vortex  line  but  may  have 
different  constants  for  different  streamlines  and  vortex  lines. 

For  an  irrotational  flow  ft  =  0  and  there  is  a  velocity  poten¬ 
tial  $ : 

v  *  V<t*  ,  (69) 

Eq.  (68)  is  valid  along  every  curve  in  the  flow  field  since 
in  Eq,  (66)  VX  *  0  everywhere.  Furthermore,  from  Eqs.  (62) 
and  (69)  we  obtain  the  Bernoulli's  equation  for  a  potential 
flow 

J*~p'  *  §■  ^  +  'ft  *  a  cons^an^  f  {IQ'S 

which  holds  throughout  the  flow  field.  The  integral  in  Eq. (70) 
must  be  computed  with  isentropic  pressure-density  relation 
because  as  discussed  in  the  previous  section,  irrotational 
motion  can  generally  be  maintained  only  by  constant  specific 
entropy  throughout  the  field. 

For  a  perfect  gas  with  constant  specific  heats,  the  equation 
of  state  can  be  written  in  terms  of  S,  p,  p  as  [Ref.  18, 

Eq.  (7-16)1 


p  «  const  pYeS^Cy 


(71) 
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Then  Eq.  (70)  becomes  [Ref.  18/  Eq.  (8-3)]  for  ^=0  and 
S  =  const  (isentropic) 

H  +  Iv2  +  rTp’lt  +  H*a  00asta»t  •  (72) 

Therefore  ir rotational  flows  are  not  necessarily  isoenergetic 
in  the  sense  that  the  total  energy  H  is  a  constant  throughout 
the  field.  Irrotational  flows  are  isoenergetic  only  if  they 
are  steady,  i.e.,  3$/3t  «  0. 

On  the  other  hand  for  a  rotational  flow  Eqs.  (67)  and  (68)  show 
that  when  the  motion  is  steady  (i.e.,  3y/3t  *  0)  the  total 

energy  H  is  a  constant  along  any  streamline  or  vortex  line 
[compare  with  the  energy  equation  Eq.  (37)].  The  difference 
is  that  for  steady  (rotational)  isoenergetic  flow  the  motion 
is  not  necessarily  isentropic,  while  for  steady  irrotational 
(isoenergetic)  flow  the  motion  must  be  isentropic.  An  example 
for  the  isoenergetic  but  nonisentropic  flow  is  the  problem  of 
steady  supersonic  flow  over  a  body  with  curved  detached  shock. 
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SECTION  4 


TRANSFORMATION  OF  CONSERVATION  EQUATIONS  IN 
TERMS  OF  ORTHOGONAL  CURVILINEAR  COORDINATES 

The  vectorial  forms  of  the  mass,  momentum  and  energy  conserva¬ 
tion  equations  Eqs.  (32),  (33),  and  (34)  can  be  expressed  in 
a  general  curvilinear  coordinate  system  by  substituting  the 
expressions  for  the  gradient,  divergence,  and  curl  operators 
in  such  a  system  given,  for  example,  in  Reference  20,  Appendix  C. 
For  practically  useful  orthogonal  curvilinear  coordinates  x^, 
x2,  x3,  with  local  unit  vectors  e^,  e2,  e3  in  the  directions 
of  increase  of  x^,  x2,  x3  the  specific  expressions  can  be 
found  in  Ref.  21,  Appendix  2. 

An  elementary  displacement  can  be  written  as 

ds  =  hjdXje^  +  h2dx2e2  +  h3dx3© 3  ,  (73) 

where  hp  h2  and  h3  are  the  scale  factors.  The  velocity  com¬ 
ponents  are  v1#  v2  and  v3  such  that 

v  -  v^  +  v^  +  v^  .  (74) 


Then  Eqs.  (32),  (33)  and  (34)  can  be  written  as  follows 
(Ref.  16,  Section  1.2.2)  : 


mass: 


momentum: 


If  +  iqKpq  (kj  hih2h3pvj)  *  0  • 

A  <pvl)  +  EEJ  (kJ  hlh2h3,lj)  + 

1  /,  shl  ,  9h2\ 

+  h^EJ  p.2  '  *22 13^7 

+  Kft  (J13  ^  '  *33  ■  £-Sl  ' 
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and  two  morons  by  cyclic  permutation 

z&zv  /t  ?t-  jsj  Ws  [<pE+p>  vj 

T^-k^^ll  =  -v-  +  G 


where  the  su  convention  has  been  used  and  9\  ^  are  the 
components  omsor  J  : 


pw  -  a  -  pw  +  pi  -  t 

S3  — —  *  «M>  MM 


They  can  be  ed  as 


with 


where 


i  *  pvivj  +  pSij  -  Tij 


»  \{V*v)6i^  +  2y  Di.  , 


5ij  @  0 


(i  ®  3) 
+  3) 


is  the  Kronecibol. 


In  Eq.  (80) 


iqtgq  hih2h3vi)  . 

l/9vl  „  V2  9hl  ,  v3  9hl\ 

TO  -EJ  Fq  +  EJ  THjj  ' 
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with  D22  and  following  by  cyclic  permutation  and 


Dij(i*j)  =  3 


*»»  i  -*■  ueLLXLi  wv-  c.  *»-•■ 


i  3vi  l  3vi  l  /  3hi 
lixj  +  hT  "  EJTT  \vi  53q  +  vj 
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.  (84 > 


SECTION  5 


CONSE  EQUATIONS  IN  CARTESIAN  COORDINATES 


To  the  Caicoordinates 


=  x,  X2  =  y,  x3  = 


(85) 


there  corrthe  scale  factors 


h^  — *  h3  "  h^  —  1  « 


(86) 


The  cons^iequations  Eqs.  (75),  (76)  and  (77)  become 
when  Eq .  (used : 


mass:  —  (pv.,)  =  0  # 

x.  j 


momentum 


5  i}  +  "5x7  (pvivj}  *  If.  ~ 


(87) 


3p  Wf  =  fi  '  (88) 


energy:  )  + 


(pE+p)Vj  - 


Tijvi  * k  Ir: 


fivi  +  0  • 

(89) 


where#  sincEqs.  (82)  to  (84) 

n  3VH 

V  •  V 


and 


n  .  i/9vi  *  !li\ 

dm  “  +  557I 


we  have #f i(80) 


(90) 


(91) 


.  .  3vi , ,  /3vi 

■ x  4ij  ^ +  “(nj 


3v,  3v 
+ 


(92) 


These  resiae  with  those  usually  found  in  the  literature 
expressed asian  tensor  notations  (see,  for  example, 

Ref.  9,  Ec?9) ) . 
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SECTION  6 


CONSERVATION  EQUATIONS  IN  CYLINDRICAL  COORDINATES 


To  the  cylindrical  coordinates 


x,  =  r,  x,  =  6,  x,  =  z 


there  correspond  the  scale  factors  (Reference  22,  Appendix  4) 


\  •  1,  h2  =  r,  h3  =  1  . 


The  conservation  equations  Eqs.  (75),  (76)  and  (77)  become 


mass: 


If  +  F  3?  (rpvr}  +  ?  JS  (pv0)  +  Jz  <PV  *  °'  (95) 


momentum:  ^  (pvr)  +  ~  ^  (ryllJ  +  f  JF  (912)  +  Jz  (913} 


-  —  3-  _  a  f 

r  22  r  ' 


Jt  (pv9)  +  r  Jr  (rff21)  +  I  W  (522)  +  Jz  {y23) 


+  r  *21  “  f0  ' 


Jt  <PV  *  I  Jr  <r*31>  +  I  W  (332> 


va33 


<3„)  -  f. 


enercj^:  ^  (pE)  +  Jr|(pE+p)vr  - 


(Tllvr+T21v0‘4'T31V2) 
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+  |  rS  [  tPE+p)Vg  -  <t12VT22V6+T32Vz> 

k  3T] 

r  70  J 

+  ^  [(PE+P)VZ  -  (T13Vr+T23V6+T33Vz). 

-  k  SI 

K  3zj 


*  frvr  +  Ve+  V2  +  G 


(99) 


where  3^  ia  given  by  Eq.  (79)  with  given  by  Eq.  (80) 
in  which, according  to  Eqs.  (82)  to  (84), 


V*v 


1  3 
r  ar 


(rvr) 


3ve  + 

ir  + 


(100) 


and 


3v  ,  /  3vq  \ 

D11  “  7F  '  °22  "  r  \TT  *  vr/  '  °33 

„  n  l  /8ve  ve  .  l  8vr) 

D12  *  d21  *  7  \-5F  •  T  +  r  T57  ’ 

„  „  i  (i  8vz  .  8ve\ 

d23  “  d32  "  7  \r  "ST  +  T?7  ' 

-  / 3v  3v  \ 

°31 53  °13  -  Ibr  +  tf]  • 


y  uoi) 
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SECTION  7 


CONSERVATION  EQUATIONS  IN  SPHERICAL  COORDINATES 


To  the  spherical  coordinates 


x,  =  r,  x,  *  6,  x, 


(102) 


there  correspond  the  scale  factors  (Reference  21,  Appendix  2) 


h^  =  1,  hj  =  r  ,  hj  =  r  sin  8 


(103) 


The  conservation  equations  Eqs.  (75),  (76)  and  (77)  become 


mass: 


H  +  ^2~^  [k  (r2sinepvr)  +  ^  (rsin0pv0) 


+  <rpv„)]  *  0  , 


(104) 


momentum: 


3  ■ 


7E  (PV  +  7r  <r2®in08ii)  +  ^  <rsine9ri2) 


+  3?  (r^i3J 


r  *22 


-  y..  m  f  , 
r  *33  r 


(105) 


A  {PV  +  [ife  (r2sin0*21>  t  (r*in6f22) 

+  (r523)]  “  ^^“^33  ”  r  *21  *  f6  '  (106) 

■k  (|>V  +  7^7  (rJsineff31)  +  &  trsineJ32) 

+  -k  «r33>]  +  5  *31  +  S  *32  -  V-  (107> 
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-i.->  \-v k-  *  .•> v> v»  k»  ->  wry*  .> i>  •»>  ■»>  \>  ■>  .-•» y»  y*  >'s y 


energy:  jr  (PE)  +  -j— —  [r2sin8  (pE+p)vr  -  ^Tnvr+T21vef T3l’V 

r  sinfi  v 

-  *  §f]j  +  TO{r3iri0(,pE+p)ve "  (T«vr+t22Te+T32V 

-  I  H])+  h  (rl(PE+P>  \  -  ‘'uVW'ssV 

•  Tzhrt  w]  =  frvr  +  f6v9  +  +  G  *  u08) 

where  7^  is  given  by  Eg.  (79)  with  given  by  Eq.  (80), 
in  which  according  to  Eqs.  (82)  to  (84), 

V*v  ■  (r2sin8vr)  +  jg  (rsin0vQ)  +  ^  (rv^)l  ,  (109) 

r*sin0  1 

and 

DU  *  "5F  '  °22  “  |(tT  +  vrj'  °33  “  rsIHF 

d12  "  d21  "  I  (tF  '  T  +  r  -wj  • 

3v-  3vq  \ 

TT  *  T*  *  cos9v^  ' 

1  / 8vr  3vtf  \ 

d31  -  d13  *  TSSSSS  +  rainQ  ‘  3in0V)  ‘ 


d23  “  °32  "  2rsOT  8in0 


hss  (B 
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SECTION  8 


COMPACT  ONE-DIMENSIONAL  FORMULATION 

The  conservation  equations  in  Sections  5,  6  and  7  for  three- 
dimensional  Cartesian,  cylindrical  and  spherical  problems 
respectively,  can  be  readily  reduced  to  those  for  various  one- 
and  two-dimensional  problems  by  setting  appropriate  velocity 
components,  force  components  and  derivatives  to  zero.  In 
particular,  if  we  consider  the  one-dimensional  problems  in 
which 


r,  Y 


-lvr(r) ' 


0, 


3(  ) 


8(  ) 
■3*T 


o,  f2  =  f3  =  o  , 

dll) 


the  conservation  equations  in  Cartesian,  cylindrical  and 
spherical  coordinates  can  be  compactly  expressed  in  one  single 
set  of  equations  (with  j  »  0,  1,  2  representing  plane,  line 
and  point  symmetry  respectively) : 


mass; 


P P,  X  >  ^l | .  j 

3t  JT 3r 


(r3pvr)  -O', 


(112) 


momentum s 


A  (pvr1  +  “J  Tr  (rJpvr2) 


.  £  -  |£  +  J- 

r  3r  3r 


3r 


-  ¥  vr  li  '  <113> 


energy:  ~  (pE)  +  -3  ^  [rjvr(pE+p)j 


£  V  + 
r  r 


1  3 

P  5? 


r3vr  j  7J  (r5vr)  +  2u  TT 


3v_ 


+  G 


(114) 
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SECTION  9 


CONCLUSIONS 


Following  a  survey  of  the  basic  equations,  the  mass,  momentum 
and  energy  conservation  equations  and  the  components  of  the 
deformation  and  stress  tensors  for  a  nonsteady  flow  of  a  com¬ 
pressible  viscous  single-phase  fluid  have  been  expressed  in 
vector  and  dyadic  notations,  transformed  in  terms  of  ortho¬ 
gonal  curvilinear  coordinates  and  specialized  for  the  three- 
dimensional  cases  in  Cartesian,  cylindrical  and  spherical 
coordinates.  These  equations  can  then  be  reduced  to  corres¬ 
ponding  equations  for  one-  and  two-dimensional  flows  in  various 
coordinates.  They  are  in  a  format  readily  applicable  to 
practical  problems. 

If  desired  the  procedure  can  be  extended  to  general  nonortho- 
gonal  curvilinear  coordinate  systems. 

This  report  also  contains  a  generalization  of  the  conservation 
equations  in  Integral  form  based  on  an  extension  of  the  Reynolds 
transport  theorem. 
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